Hopping Models of Charge Transfer in a Complex Environment: New Class of 
Coupled Memory Continous-Time Random Walks 
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Charge transport processes in disordered complex media are accompanied by anomalously slow 
relaxation for which usually a broad distribution of relaxation times is adopted. To account for 
those properties of the environment, a standard kinetic approach in description of the system is 
addressed either in the framework of continuous-time random walks (CTRW) or fractional diffusion. 
In this paper the power of the CTRW approach is illustrated by use of the probabilistic formalism 
and limit theorems that allow to predict the limiting distributions of the paths traversed by charges 
and to derive effective relaxation properties of the entire system of interest. Application of the 
method is discussed for non-exponential electron-transfer processes controlled by dynamics of the 
surrounding medium. 



PACS numbers: 05.40-|-j, 82.20.Fd, 87.10. -)-e 



I. INTRODUCTION 

The stochastic formulation of transport phenomena in 
terms of a random walk process, as well as the descrip- 
tion via the deterministic diffusion equation are two fun- 
damental concepts in the theory of diffusion in complex 
systems. The best known examples are charge trans- 
port in amorphous semiconductors, rebinding kinetics in 
proteins, polarization fluctuations in inhomogeneous sol- 
vents, diffusion of contaminants in complex geological 
formations and diffusion of pollutants in large ecosys- 
tems. In all realms mentioned above, the complex struc- 
tures, characterized by a large diversity of elementary 
units and strong interaction between them, exhibit a non- 
predictable or anomalous temporal evolution. The possi- 
bility of the dual description of the anomalous dynamical 
properties of such systems, based either on the random 
motion or on the differential equations for the probability 
density functions, has been considered in literature since 
the late 60s and gave rise to an extensive list of developed 
models ^2, 3]. 

In this paper we demonstrate the power of the math- 
ematical tools underlying the concept of a continuous- 
time random walk (CTRW) by showing how the tool can 
be generalized to handle complicated situations such as 
diffusion-reaction schemes in complex system. The no- 
tion of the CTRW, a walk with a waiting time distri- 
bution governing the time interval between subsequent 
jumps of a random walker, has been introduced by Mon- 
troU and Weiss Q. The distribution of waiting times 
may stem from possible obstacles and traps that delay 
the particle's motion and in consequence, introduce the 
memory effect into the kinetics. Especially fascinating 
in this approach was the idea of an infinite mean time 



between the jumps as in such characteristic time 

scale of the process looses its common sense. This novel 
concept has been used by MontroU and Scher to give 
a first explanation of experiments measuring transient 
electrical current in amorphous semiconductors. Since 
then the CTRW formalism has been successfully applied 
to describe fully developed turbulence, transport in frac- 
tal media, intermittent chaotic systems and relaxation 
phenomena. The common feature of the above men- 
tioned applications is that they exhibit anomalous dif- 
fusion manifested by a non-Gaussian asymptotic distri- 
bution (propagator, diffusion front) of a distance reached 
at large times. 

At the level of the CTRW modeling, the diverging mean 
waiting time leads to a subdiffusive motion with the mean 
square displacement growing as < r^(t) >(x with 
< a < 1. When applied to the theory of Brownian 
motion, the CTRW scenario leads to the fractional dif- 
fusion equation 's', that can be treated on an equal 
footing with the framework used for systems with nor- 
mal diffusion. 

Usually, in applications of the CTRW ideology, the anal- 
ysis of the asymptotic distribution is presented within 
the approach that is based on a formal expression for 
the Fourier-Laplace transform of the propagator, or oth- 
erwise, use of the fractional calculus is required 01 as 
a legitimate tool. Here, we present an approach to a 
random walk analysis which is based directly on the def- 
inition of the cumulative stochastic process. Our aim is 
to show that despite of the extensive studies on CTRWs 
and their long history in physics, the powerful tool of the 
limit theorems hidden behind the derivation of limit- 
ing distributions, has not been fully explored yet. We em- 
phasize the possibilities of applications of that scenario 
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in stochastic modeling of physical systems, in particu- 
lar, in description of the charge transport in disordered 
materials. 

A starting point in the CTRW analysis is the defini- 
tion of a total path R{t) of a particle traversed up to 
the time t in accumulating number L{t) of jumps of a 
length Rj. The number of jumps exerted in (generally 
random) time t can be defined either directly by assum- 
ing a specific counting process L{t) (with e.g. Poisson, 
negative binomial, geometric, etc. count distribution) 
or indirectly - by assuming the distribution of waiting 
times Tj between the jumps. In both approaches, un- 
der certain assumptions concerning the distribution of 
jumps Rj and the distribution of waiting times Tj (or 
number of jumps L{t)), the asymptotic distribution of 
the total path R{t) reached up to time t can be obtained 
by applying limit theorems of the probability theory. In 
contrast to the more popular Tauberian analysis of the 
Fourier-Laplace transform of -R(i), such an approach pre- 
cisely identifies classes of possible limiting distributions 
and offers an easy-to-follow scheme of generating various 
limiting results. 

The dual description of the random walk and the rela- 
tionship between the results obtained in both cases are 
exemplified in this work by discussing the biologically 
relevant charge transport processes in disordered media. 

The paper is organized as follows: We begin in Section 
II with a brief discussion of models of non-exponential di- 
electric relaxation and their relation to solvent (medium) 
dynamics infiuencing the rates of the long-range electron 
transfer. Further, as a generalization of the McConnell 
formula we incorporate medium fluctuations in the ex- 
pression for the electronic transfer matrix. Its form is 
analyzed in terms of an exponential of a sum of indepen- 
dent and identically distributed (i.i.d) random variables 
with a random number L of virtual jumps between the 
donor and acceptor sites. By assuming the deviations 
from equilibrium of the atomic coordinates of a given 
pathway to be random contributions to the sum, we are 
able to investigate asymptotic forms of the tunneling ma- 
trix elements. Section III poses the problem in terms of a 
standard CTRW scenario which is generalized (Sections 
IV and V) for random walks subordinated to a compound 
step-counting process. Main results and conclusions of 
the analysis are presented in Sections VI and VII. 



II. CHARGE TRANSPORT IN A COMPLEX 
ENVIRONMENT 

Charge transport processes determine a variety of phe- 
nomena in physics, chemistry and biology. The study of 
the phenomenon has gradually developed together with 
general progress in theoretical physics and in fast high 
resolution spectroscopy, so that contemporary research 
deals nowadays with a broad class of systems, materi- 
als and environmental conditions. Of particular inter- 
est are the processes taking place in disordered mate- 



rials [3, such as amorphous semiconductors, randomly 
arranged molecular wires, glasses or biological proteins 
where the charge transfer processes form the elementary 
steps in energy transport and production in almost all 
living cells. In all those cases, the actual transport pro- 
cess is coupled to local polarization fluctuations of the 
environment. For the situations that the relaxation of 
the polarization fluctuations of the surrounding medium 
has a simple "close-to-equilibrium" exponentially decay- 
ing form, the main energetic contributions to the charge 
transfer process come from the reorganization energy of 
the medium [ol Il0||. In contrast, many observed charge 
transport processes, like electron transfer (FT) in com- 
plex solvents and proteins 
or gating kinetics of biological channels |20j , 
exhibit non-exponential kinetics resulting from the com- 
plex response to the interfering medium. A classical ex- 
ample are higher alcohols, for which the frequency de- 
pendent dielectric permittivity takes on a Cole-Davidson 
(CD) [13 form: 
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(1 + iujTp)'r 



dt 



(2.1) 



with < 7 < 1; (j){t) standing for the correlation function 
of polarization fluctuations and Tp indicating a reciprocal 
proportional to the peak frequency of the dielectric loss. 
In the electron transfer (FT) theory, the time-correlation 
function (j){t) is related 0, |^ to the Coulombic po- 
tential energy difference for a given configuration of all 
solvent (intervening medium) molecules in the states of 
reactants and products: 



0(i) = {{Sl^Ef) ^ {SAESAE{t)) 



(2.2) 



with AE{t) identified with a complex dynamic "reac- 
tion coordinate" describing the transfer. In a convenient 
dipole-approximation for medium molecules, the poten- 
tial energy difference AE would be given by 

^E = -J drP{r)[Epir-rp)-ER{r-rn)] (2.3) 

where P(r) stands for the medium orientational polariza- 
tion at position r. For a solvent in which the dipoles of 
the dielectric medium relax with a single relaxation time 
Tp, the complex dielectric permittivity Fa. H2.1|l is given 
by the Debye (D) function 
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(2.4) 



with (j){t) expressed in terms of a single exponential func- 
tion with a decay time Tp. Other, equally likely fitted 
expressions [Tsl Il6j exploited in dielectric spectroscopy 
of polymers and disordered solids estimate relaxation of 
(t){t) by use of the Cole-Cole (CC) [H formula 
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or the Havriliak-Negami (HN) [T^ function: 
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(2.6) 



where < a < 1 and < 7 < 1 are parameters determin- 
ing the characteristics of the dielectric relaxation with a 
representing the width and 7 the skewness of the dis- 
tribution of relaxation times [iSj. Although the generic 
physical reasons for anomalous relaxation in complex sys- 
tems are still under debate, both - static models based 
on the inhomogeneity of the medium, - as well as the dy- 
namic models, describing complex local dynamical pro- 
cesses have been successfully employed to describe relax- 
ation behavior of fluctuations in such systems. In partic- 
ular, the studies on the effect of protein dynamics on bio- 
logical ET 0,|il|2l|2l|2| have demonstrated sensitiv- 
ity of the long distance tunneling mediated by the protein 
matrix on atomic configurations of the surroundings and 
pointed out possibility of an electron of emitting or ab- 
sorbing phonons from the medium that would effectively 
result in an inelastic ET processes. In numerous chemical 
and biological examples of the ET reaction ^25, '23'| , a 
single electron is tunneling in an inhomogeneous medium 
over large distances of several angstroms. The interven- 
ing medium can be either a protein backbone or a se- 
quence of cofactors embedded in a protein matrix. Due 
to a large separation between the donor and acceptor, 
direct electronic coupling between the chromophores is 
negligible, rendering thus the question on the effect of 
medium on enhancement of the electronic coupling [17| . 
A possible realization of the long-distance ET process is 
a transfer mediated through the medium which acts as 
a bridge providing virtual states for the tunneling elec- 
tron [25| . Within the nonadiabatic-reaction scenario cor- 
responding to a weak electronic coupling Tda between 
the state of reactants D and products A, the expression 
for the rate reads 



27r 

kET = ^TU{FC) 



(2.7) 



where (EC) is the Franck-Condon nuclear factor associ- 
ated with the nuclear modes activation barrier. In a con- 
ventional theory the Condon approximation is assumed, 
i.e. the electronic coupling Tda is viewed as independent 
of the coordinates of the medium. To account for thermal 
fluctuations of the bridge or random intervening medium, 
the electronic coupling has to be a function of the modes 
of the medium. The simplest expression that can be pro- 
posed in such a case is the Mc Connell formula 0, 
which for a case of a linear bridge consisting of L orbitals 
leads to the tunneling matrix Tda 



DA 



n 



(2.8) 



with e — ej being the energy difference between the tun- 
neling energy and the energy of the bridging orbital j, 



L standing for the number of virtual jumps performed 
along the path and f3ij denoting couplings between di- 
rectly overlapping atomic orbitals of neighbouring atoms 
within the bridge. The above formula constitutes the 
essential part of the ET pathways models JJ^ |2J| in pro- 
teins, where the calculation of the effective electronic cou- 
pling is based on a general assumption that the electron 
wave function decay is softer for propagating through 
a chemical bond than through space jump. Since the 
coupling coefficients /3 are exponentially decaying func- 
tion of the distance between subsequent medium centers 
(atoms), the effective tunneling matrix can be recast in 
the form 




(2.9) 



where rj are fluctuations of the atomic coordinates of 
the bridge, aj are constants characterizing strength of 
the coupling to a particular bridge mode j and cor- 
responds to the average, equilibrium tunneling matrix. 
Such a representation of the effective tunneling matrix 
allows to use the notion of the continuous time random 
walk (CTRW) as a very convenient mathematical tool to 
analyze the decay with time of the donor state occupa- 
tion density. In fact, the latter is commonly described by 
systems of phenomenological balance equations 



dt I P2{t) 
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P2{t) 



(2.10) 



which relate the decay of the donor (acceptor) popula- 
tions to the state-relaxation rate constants Note, 
that corresponding populations -Pi, 2(0 ^.ny of the elec- 
tronic states (reactants or products) are dynamical quan- 
tities usually measured in the electron transfer kinetic ex- 
periment and are obtained by integrating the polarization 
energy dependent populations p{E^ t) over configuration 
variable E{t): 
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(2.11) 



In a standard ET theory approach |21|, |2^ |2J| after as- 
suming a disentanglement of reactive tunneling from the 
dynamics of diffusion, the elements of the evolution ma- 
trix Eg. H2. 10(1 have the form of 



k+ 



1 + ^na/^d + ^na/^d 



(2.12) 



where k^A describes the crossing (nonadiabatic) kinet- 
ics and ku is the rate constant of the diffusion in the 
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reactants' (products') basins. However, in a more gen- 
eral situations, where the matrix entries in Ea. H2.10|) 
are represented by time dependent functions, the re- 
distribution of populations and consequently, the relax- 
ation of clectron-donoring (accepting) states may follow 
a non-exponential law. Accordingly, the frequency char- 
acteristics of dielectric susceptibility x(a;) connected to 
the temporal relaxation function of the induced state- 
polarization P{t) = eox(w)E(t) = eo{e* — l)E(t) where 
E(t) = Ege""^* and the functional character of the di- 
electric permittivity e*{uj) may be inferred from the anal- 
ysis of relaxation of state populations in a frequency do- 
main 

poo 

X{u) = / e-*'^*d(-Pi.2(t)) (2.13) 

"'0 

In the forthcoming sections we present a dynamic frame- 
work which, within the CTRW cenario, leads to the em- 
pirically observed non-exponential relaxation dynamics. 



and 

cT{t) = L,/2{t; 1, 1,0) = ^t-3/2e-i/2t^ ^ > q, (3.4) 
V 27r 

i.e. the model describes a charge moving only in one 
direction with a Poisson number of jumps and the Levy- 
Smirnov distribution Li/2it; 1,1,0) of waiting times. The 
probability density function p2 in Eq. (|3.2|l is then given 
by a convolution of independent probability densities 
(El 

P2(r) = Li/2(r;j2,l,0) (3.5) 
and leads to a marginal probability density 

.(T).^e-^f:i^,e-^ (3.6) 



III. PATHWAY ANALYSIS OF ET REACTIONS 

As discussed above, with the matrix elements of a par- 
ticular path Td^, the nature of disorder may be analyzed 
in terms of fluctuations in couplings or, alternatively, in 
contributions Rj to the total distance R traversed by a 
charge: 

R = ^ Rj>n (3.1) 
j=o 

Here L{TZ) stands for a random counting process describ- 
ing a (random) number of forward steps exerted by a 
particle before reaching a distance TZ. Note, that such 
a formulation is identical with the assumption of a one- 
dimensional, biased (directed) random walk performed 
in an amorphous medium under the influence of a strong 
external field. The time T that particle needs to reach a 
fixed distance TZ is given by T = '^fl'^^ Ti with Ti being a 
time spent by a hopping charge at the location ^ Rj . 
By means of a conditional probability, the probability 
density p{T)fi for the distribution of times T to reach a 
distance TZ reads 

p{T)n = 

oo /L(n) \ 

= ^Pl(i(7^) = j)p2 E = ^l^(^) = (3.2) 

where p2 stands for the probability distribution that the 
elapsed time is T provided exactly j steps have been per- 
formed to reach distance TZ and pi is probability to make 
j steps over the distance TZ. Here we assume, that the 
length of a given jump, as well as the waiting times elaps- 
ing between two successive jumps are drawn as indepen- 
dent random variables with densities 

p{r) = /3e-'3^ r > 0, (3.3) 
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FIG. 1: Density distribution function p{R)t of distances tra- 
versed up to the time T in a one-dimensional random with 
T{t) given by a Levy-Smirnov distribution Eg. 13.41 
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FIG. 2; Density distribution function of times T needed to 
reach a distance 7?. in a one-dimensional random walk with 
a Levy-Smirnov distribution Eq.((23J of waiting times for a 
jump. 

Quite similarly, one can calculate p{R)t with the con- 
dition T — X^jlT^"''^ -^J — ^'^^ P'^ given by a gamma 
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density distribution resulting from a sum of j indepen- 
dent and exponentially distributed random variables: 



PiRh 



E 

n=0 



m)\-0n 2 



n+l 



e 2T (fy 



(3.7) 



Thus the sketched problem describes a 1-dim diffusion 
among traps with a broad (asymptotically, heavy-tailed) 
distribution (j[t) of trapping times, see Figures 1 and 2. 
This type of the "annealed" CTRW has been extensively 
studied in literature Q, S 111 IH A 3-dim ana- 

logue with a broken unidirectionality of the transfer cor- 
responds to a weak external field approximation. Fig- 
ure |31 displays results of a 3-dim CTRW computer simu- 
lation: direction of a jump has been generated by sam- 
pling spherical coordinates 9 and (p from uniform distri- 
butions defined on intervals (0,^) and (0, 27r), respec- 
tively. Quite arbitrarily, the positive direction of the z- 
axis have been favoured by sampling the coordinates of 
a point ((/), 0, —z) by using those for (0, 0, z) and switch- 
ing from z to —z with a probability 3/10. The upper 
panel of Figure |21 refers to the distribution of times T 
required to reach the distance R when the distribution 
of waiting times for the subsequent jump was assumed in 
a form of a skewed (heavy-tailed) Levy-Smirnov distri- 
bution describing higher probability of long trapping. In 
contrast, the lower panel represents results of simulations 
for a 3-dim CTRW with a preferential short mean time 
of waiting (a WeibuU distribution) for a particle release 
from the trap: 



a{t) = ai" 



(3.8) 



where < a < 1. In a forthcoming section we investigate 
an extension of such an approach for the family of com- 
plex CTRW incorporating effects of clustering of carrier 
trapping sites in a disordered material under the study. 



IV. CTRW SUBORDINATED TO A 
COMPOUND COUNTING PROCESS 

As briefly mentioned in the preceding sections, the 
CTRW generalizes a simple random walk by implement- 
ing a random waiting time between jumps ^1„ JJ. This 
stochastic process is defined by the total distance R{t) 
reached by the particle at time r > if the movement is 
generated by a sequence {{Rj,Tj), j — 1,2, . . .} of jump 
parameters with Rj specifying both the length and the 
direction of the jth jump and Tj > denoting the wait- 
ing time for the next jump. We assume moreover that 
the sequence {Rj, Tj ),_;/ = 1, 2, ... is formed by i.i.d. ran- 
dom vectors. The cumulative distance R(t) may be 




FIG. 3: Density distribution function of times T needed to 
reach a distance 7?, in a three-dimensional CTRW with a Levy- 
Smirnov Eg. 13.41 distribution (upper panel) and a Weibull 
distribution Ea. l|3.8l ('lower panel) of waiting times for a jump. 



(cf. Eq. ^'S.l\\ ) expressed as 



where 



r{t) = j2rs 



/(t) = min{A: : ^ > r} 



The renewal process 




1,2,.. 



(4.1) 



(4.2) 



(4.3) 



represents the instants of time at which subsequent jumps 
occur; the process {v{t),t > 0} counts the jumps. 
Hence, the introduced CTRW {R{t),t > 0} is the 

[r] 

discrete-time random walk {'^ Rj,T > 0} subordinated 

to the renewal counting process {t/(r),r > 0}. It is 
decoupled if random variables Tj and Rj are indepen- 
dent; otherwise it incorporates statistical dependence be- 
tween time and space steps. In general, the properties of 
{R{t),t > 0} are unknown. Nevertheless, the behav- 
ior of the CTRW for long time, or, equivalently, of the 
rescaled process R{t / St) / f {6t) , (where fiSr) is an ap- 
propriately chosen rescaling function) for the characteris- 
tic time scale St decreasing to 0, can be determined quite 
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well. Systematic studies of the limiting total-distance dis- 
tributions for the one-dimensional CTRWs have pointed 
on few possible distributions . Introducing a new 

class of coupled memory CTRWs with random-sum form 
of time/space steps, much broader class of possible limit- 
ing distributions can be obtained. Namely, let us consider 
the case when 



Mi 



Mi 



Rj = ^ ARjk , =^ ^Tjk, 



(4.4) 



fc=i 



fe=i 



where sequences {ARjk,j,k 



1,2,...} and 



{ATj).,j,k — 1,2,...} of space/time random spans 
are independent and each of them consists of i.i.d. pos- 
itive random variables, and {Mj,j = 1,2,...} is a 
sequence of i.i.d. positive integer- valued random vari- 
ables. The nondegenerate distribution of Mj, the 
number of summands, (i.e. the case when random vari- 
able Mj takes at least two different values with positive 
probabilities), provides a stochastic dependence between 
time and space steps. As a consequence, the CTRW 
resulting from the space/time-stcp family Ea. H4.4|l is 
usually coupled. Moreover, it has an equivalent random- 
sum form since the total distance R{t) has the same 
distribution as the following random sum of distance 
increments 



(4.5) 



k=l 



where the random index L{t) has the form 

L{t)^ J2 M,, (4.6) 
with the random number of summands defined as 

n 

N{K{t)) = min{n : > K{t)} 



for 



K{t) = min{n : ^ ATik > t}. 



(4.7) 



(4.8) 



fc=i 



Note, that the above definition of the CTRW process (cf. 
Ea. (|4.5|) ') does not assume exphcitly any precise relation 
between steps Ri and time lapses T^. Instead, possible 
coupling between those two random variables is incorpo- 
rated via the counting process Eq. (|4.6I) summing contri- 
butions to the cumulative distance R{t) traversed during 
the walk. This particular construction of the random sum 
Eq. (|4.5() allows calculating statistical properties of their 
distribution by use of the "nested" moment-generating 
function $(s) =< s^'-'^^ > which in this case fulfills the 
equation 

^Ris) = <I>L(<i'Afli, (s)) = $ Jv(1>A/(*Ai?,, is))) (4.9) 

where the last equation follows by applying the technique 
of conditional averages. 



V. DIFFUSION FRONT. LIMITING 
DISTRIBUTIONS 

Let us consider now the rescaled total distance 
RM{t/ST) 



RSrit) 



fiSr) 



(5.1) 



where St is the characteristic time scale, and J{5t) is a 
rescaling function chosen appropriately. Applying limit 
theorems '30] one can evaluate the limiting position R{t) 
of the rescaled total distance reached as 5t 0. The 
characteristics of R{t) depends on assumptions set on 
the distributions of the variables ARjj,, ATj^) and Mj. 
Below, following the regularization scheme Ea. l5.1|) along 
with a more detailed presentation and proofs included in 
ref.fs^,!^ , we discuss briefly some examples which might 
be of practical use in modeling relaxation phenomena in 
disordered materials. We consider the case when ARjk 
and ATjk are independent and both positive. 

(a) Let us first assume that both ARjk and ATjk have 
heavy-tailed distributions with c = ci and c ~ C2, 
respectively, and the same r — a. (We say that the 
distribution of a positive random variable X has a 
heavy tail if for some c > and < r < 1 



lim 



Pv{X>x) 



1. 



(5.2) 



In consequence, the expected value {X) is infinite.) 

• If the distribution of the random number Mj 
has a heavy tail with some c > and r = 7, 
then for any i > 

1/q 

. (5.3) 



R{t/5T) _d 
1/5t 

Here " 



d 

reads "tends in distribution" , and 



"=" denotes the equal distributions. 
• If the numbers Mj have a finite expected value 
{{Mj) < 00), then for any t > 



l/Sr sr^o"^^^' - A So,- 



(5.4) 



(b) Assume that the expected values of both ARjk and 
ATjk are finite, and {ARjk) = ci, {ATjk) = C2. 

• If the distribution of Mj has a heavy tail with 
some c > and r = 7, then for any t > 



R{t/ST) 
I /St 



(5.5) 



If (Mj) < 00, then for any t > 



R{t/ST 
1/St 



a.s 



m 
1 

A 



(5.6) 



Here ' 



*" reads "tends with probability 1". 
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The random variables B^, Sa, and in 15. 3() - (|5.6() 
are as follows: 



• Sa and are identically distributed according to 
the completely asymmetric a-stable law such that 

• is distributed according to the generalized arc- 
sine distribution with parameter 7 (i.e., the beta 
distribution with parameters p = 7 and (7 = 1 — 7) 
given by the density function 



r(7)r(i-7) 




(;'y-i(l-a;)^T for < x < 1 



otherwise; 



(5.7) 



• for any < 01,7 < 1 the random variables Bj, Sa, 
and Sa are independent. 

Armed with the above results, we are now in position to 
discuss properties of CTRW paths generated under men- 
tioned constraints with the application of the formalism 
in the analysis of the relaxation responses in disordered 
materials. 



VI. EMPIRICAL AND PHENOMENOLOGICAL 
RELAXATION RESPONSES 

Relaxation in amorphous solids and ET processes in 
disordered molecular media represent nowadays inten- 
sively investigated subjects both in experimental and the- 
oretical physics 0.liljllMl2l)l2l)l23- particular, a key 
probe of electron dynamics in disordered systems is the 
time of flight experiment (TOF) for the drift mobility. In 
the experiment, the thin film sample is located between 
two blocking contacts across which is maintained a po- 
tential drop, and a laser flash is used to create carriers 
that wander towards an appropriate electrode. During 
their drift through the sample, the electrons and holes en- 
counter a variety of traps that affect their motion. The 
experiments show that in the disordered materials, the 
registered transient current follows an algebraic decay 
I{t) w In contrast, for Gaussian transport pro- 

cesses, the charge carriers move at a constant velocity and 
after a transient time, depending on the thickness of the 
sample and the applied external field, they become ab- 
sorbed. In consequence, for normal transport processes 
typically observed in ordered materials, the current is 
given by a step-like function, whereas in disordered me- 
dia the current I{t) adheres to a universal (independent 
of the applied field and sample thickness) scaling curve. 
Similar conclusions are drawn from the ultrafast pump- 
probe laser spectroscopy and spectral hole burning ex- 
periments which are well advanced techniques used for 
nanostructures and comprise nowadays a standard tool 
to determine fast carrier dynamics and spectral and spa- 
tial diffusion of the carriers. A wide-ranging experimen- 
tal information resulting from the latter [isl l33 | has led 



to the concept that the classical phenomenology of re- 
laxation processes breaks down in complex materials. It 
has been found that the Debye behavior 12.4|l is hardly 
ever found in nature and that for many dielectrics the 
deviations from it may be relatively large 0, [l5| . 
For a long time a major effort has been diverted to 
a purely qualitative representation of the shape of the 
non-Debye dielectric functions in terms of certain math- 
ematical expressions without, in any way, going into a 
physical significance of these representations. As pointed 
out by experimental studies almost all dielectric data are 
characterized well enough by a few empirical functions 
[Til IT^ IT^ Hsl IT^ . The most popular analytical expres- 
sion applied to the complex susceptibility or permittivity 
data is given by the Havriliak-Negami function (|2.t)|l . For 
a = 1 and 7 < 1 , formula (|2.t)|l takes the form (|2.1|) of the 
Cole-Davidson function; for 7 = 1 and a < 1 it takes 
the form 1)2. 5|l of the Cole-Cole function, and for a = 1 
and 7=1 one obtains the classical Debye form (|2.4|l . 
Let us note that time-domain relaxation function (j){t) 
corresponding to formula H2.6|l has the following series 
representation 



, /,N 1 (-l)"r(7-|-n) ,xa(7+«) 

^^^^^'"Sr(7).!r(i + .(7 + n)) ^"^^'^ ■ 

(6.1) 

referred to the generalized Mittag-Leffler distribution. In 
case of the CC function the series representation Ht).l|l is 
simplified to 



^ r(l + a[l + n)) 



(6.2) 



corresponding to the Mittag-LefSer distribution. The CD 
relaxation function is referred to the tail function of the 
gamma distribution with the scale parameter and the 
shape parameter 7, given by the density function 



^^(^)^ JiT%Mr''e~->'* for t>0. 







otherwise. 



(6.3) 



In order to derive relevant relaxation functions resulting 
from cases 1)5.311 - (|5.6|l . considered above, we use the 
following relations; For any < a < 1 



and for any < 7 < 1 we have j31| 



(6.4) 



(6.5) 



where S\ = S-^ = B\ = \\ G\ = E \s exponentially dis- 
tributed with mean 1 ; and for 7 < 1 the random variable 
is distributed according to the gamma distribution 
with the shape parameter 7 and the scale parameter 1. 
Assuming that for any < 01,7 < 1 the random variable 
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Table 1 . Special cases of the frequency-domain 
responses (|6.8|l . 
(co is a positive constant.) 



Assumptions 


Response 


ARjk 




Mj 


heavy tail 
r = Q 


heavy tail 
r = a 

C=Co 


heavy tail 
r = 7 
OO 


Havriliak-Negami 
Q, 7 < 1 


heavy tail 

r = a 
c = ujpCo/k 


heavy tail 
r = a 

C=Co 


{M,><oo 


Cole-Cole 

a< 1, 7 = 1 


ojpCo/k < oo 


(AT.O = 
Co < oo 


heavy tail 
r = 7 
OO 


Cole-Davidson 

Q = l, 7<1 


ojpCo/k < oo 


(AT,0 = 
Co < oo 


{M,><oo 


Debye 

Q = 1, 7=1 



is independent of Sa , and by using the conditional- 
expectation tools, from formulae Ht).4|) and Ht).5|l one ob- 
tains 



-k 



1/a 



Pr(Gy"5„ > k) (6.6) 



and hence for R{t) of the form H5.3|l - H5.6|l the corre- 
sponding relaxation function (j){t) equals 



-fci?,(t)\ 



Pr((A/fc)Gy"5„ > i). (6.7) 



where k is an appropriate positive constant. For such 
a relaxation function we have the frequency domain re- 
sponse of the general form 



= /p-«('^/'^p)G^ 



LUp = k/A, 



(6.8) 



which includes the Havriliak-Negami with its special 
cases (see Table 1). 



VII. CONCLUSIONS 

Relaxation processes deviating from the usual exponen- 
tial behavior in time domain (and a classical Debye 
form in the frequency domain) occur in many physical, 
chemical and biological systems, such as supercooled liq- 
uids, viscoelastic solids, polymer melts and porous me- 
dia, membranes and liquid crystals. They are usually 



described by some mathematical functions related to the 
fractional-order differential equations. Properties of re- 
laxation processes in complex media are also important 
topics in the long-range transfer of an electron in poly- 
meric and various biological materials 0, Is^. Issl ISq. 
The kinetics of the charge transport in those media is 
determined by the nature of electronic coupling which 
for the long distances is mediated by sequential overlaps 
of atomic orbitals of the donor, the intervening medium 
(bridge), and the orbitals of the acceptor. Internal ran- 
dom motion of the medium may result in fluctuations of 
the tunneling barriers between subsequent transfer states 
and modulate the electronic coupling. Such an effect is 
formally due to the dependence of the electronic coupling 
on the nuclear coordinates of the medium and may ex- 
press possible deviation from the usual Condon approx- 
imation ll^,!!^- The influence of factors arising from 
the static and dynamic fluctuations in electronic coupling 
have been discussed in a number of papers (|24j, 18], 3^1 
and references therein) related to the electron transport 
in proteins and polymers. All of them have claimed ex- 
istence of nontrivial effective coupling resulting from av- 
eraging over environmental disorder. 
In this paper we have extended previously proposed 
model |36j which accounts for fluctuations in long dis- 
tance charge transfer hopping. Our approach is based on 
CTRW formalism in the representation of random sums 
of i.i.d elements which are deviations from equilibrium 
of the atomic coordinates of the intervening medium. In 
contrast to other works 5, 6, 27], the present approach is 
based on renewal theory and the limit theorems for ran- 
dom sums |3,l2Sl2Sl2fll23 of jumps instead of Tauberian 
theorems for the two-dimensional Laplace-Fourier trans- 
form. Dispersive kinetics, as discussed here, appears in 
situations when the environmental fluctuations become 
comparable or slower in the time-scale of the overall tran- 
sition from the reactants' to products' states. In ordinary 
solvents, the medium relaxation is complete in several pi- 
coseconds. In contrast, in highly viscous liquids the re- 
laxation modes can determine the rate of the fast ET ki- 
netics. Theoretical treatment of such processes is based 
on analysis of slow, reorientational motions of the sys- 
tem which influence survival of populations in given elec- 
tronic states. Our analysis maps the complex medium 
reorientation kinetics on a relaxation process of measur- 
able reactants' and products' populations for states be- 
fore and after the charge transfer. Diffusion of a charge 
(and response kinetics of the local polarization energy) 
can be then investigated by use of the generalized relax- 
ation functions defined as moment generating functions 
for coupled memory random sums. A non-exponential, 
non-Debye decay of relaxing modes is a basic feature 
arising in such an approach and can be traced to the 
probabilistic pattern of contributing lengths of steps and 
waiting times in the charge hopping process. 
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